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1 Introduction 

Symmetries play an important role both in physics and mathematics. They are described 
by transformations leaving structural relations unchangeable. Their importance range from 
fundamental and theoretical aspects to concrete applications, having profound implications 
in the dynamical behavior of the systems and in their basic qualitative properties. Knowledge 
of symmetries result on a deep insight about the inner structure of a system, and permits 
to apply the conservation laws to the investigation of the objects, i.e. to link the invariance 
principles with the conservation laws. This interrelation includes three classes of the most 
fundamental principles of physics: symmetry, conservation, and extremality. 
When a closed system is characterized by a quantity which remains unchangeable in the course 
of time, no matter what kind of processes take place in the system, such quantity is said to be 
a conservation law. Some fundamental conservation laws include the conservation of energy, 
impulse, momentum impulse, motion of the centre of gravity, electrical charge, and others. 
All physical laws are described in terms of differential equations (equations of motion). The 
conservation laws represent the first integrals of the equations of motion and are important 
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for three reasons. Firstly, the task of solving the equations of motion explicitly is not always 
possible and knowledge of the first integrals may considerably simplify that task. Secondly, 
often there is no more necessity to solve the equation of motion, as the useful information 
is contained in the conservation laws. Thirdly, the conservation laws have a deep physical 
meaning and can be measured directly. 

All basic differential equations of physics (i.e. the equations of motion of physical systems) 
have a variational structure. In other words, the equations of motion of a physical system 
are the Euler-Lagrange equations of a certain variational problem. It turns out that the 
conservation laws are the result of the invariance of the action with respect to a continuous 
group of transformations, given by some symmetry principle. The more general expression of 
the interrelation symmetry/variational structure/conservation, is given by Noether's theorem. 
Noether's theorem asserts that the conservation laws for a system of differential equations 
which correspond to the Euler-Lagrange equations of a certain variational problem, come 
from the invariance of the variational functional with respect to a one-parameter continuous 
group of transformations. The group of symmetry transformations requested by Noether's 
theorem depend, of course, on the physical properties of the system. We refer the interested 
reader to [7]. 

Conservative physical systems imply frictionless motion and are a simplification of the real 
dynamical world. Almost all systems contain internal damping and are subject to external 
forces. For non-conservative dynamical systems, i.e. in the presence of non-conservative 
forces (forces that do not store energy and which are not equivalent to the gradient of a 
potential), the conservations law are broken so that the standard Lagrangian or Hamiltonian 
formalism is no longer valid for describing the behavior of the system. Methodologically, 
Newtonian dissipative dynamical systems are a complement to conservative systems, because 
not only energy, but also other physical quantities as linear or angular momentums, are not 
conserved. In this case the classical Noether's theorem ceases to be valid. However, it is still 
possible to obtain a Noether-type theorem which covers both conservative (closed system) and 
nonconservative cases [H [5]. Roughly speaking, one can prove that Noether's conservation 
laws are still valid if a new term, involving the nonconservative forces, is added to the standard 
constants of motion. 

In order to better model non-conservative dynamical systems, a novel approach entitled Frac- 
tional Action-Like Variational Approach (FALVA) has been recently introduced [11 [3]. This 
approach is based on the concept of fractional integration. Fractional theory plays an impor- 
tant role in the understanding of both conservative and non-conservative behaviors of complex 
dynamical systems, and has important physical applications in various fields of science, e.g. 
physics, material sciences, chemistry, biology, scaling phenomena, etc. In [2^, '3] Riemann- 
Liouville fractional integral functionals, depending on a parameter a but not on fractional- 
order derivatives of order a, are introduced and respective fractional Euler-Lagrange type 
equations obtained. In [8], Jumarie uses the variational calculus of fractional order to derive 
an Hamilton-Jacobi equation and a Lagrangian variational approach to the optimal control of 
one-dimensional fractional dynamics with fractional cost functional. In this paper we extend 
the results of [1] to more general FALVA problems with higher-order derivatives. 

2 Preliminaries 

We begin by collecting the necessary results from [HE]. In 2005 El-Nabulsi (cf. [2]) introduced 
the FALVA problem as follows: 
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Problem 2.1. Find the stationary points of the integral functional 

i[qi-)] = -l- [\ie,q{e),q{e))it-er-'d9 (2.1) 

r(a) Ja 

under the initial condition q(a) = qa, where q = T is the Euler gamma function, < a < 1, 
9 is the intrinsic time, t is the observer time, t ^ 9, and the Lagrangian L : [a, b] xR" xM" M 
is a C'^ -function with respect to all its arguments. 

Along all the work we denote by diL, i = 1,2, 3, the partial derivative of L{-, ■, •) with respect 
to its ith. argument. 

Theorem 12.21 summarizes one of the main results of [2]. 

Theorem 2.2 ([2j). If q{-) is a solution to Problem \2.1\ (i.e., q{-) is a critical point of the 
functional (j2.ip ). then q{-) satisfy the following Euler-Lagrange equation.- 

d2L {9, q{9),q{9)) " ^^sL {9, q{9),q{t)) = ^^3^ q{e)A{e)) • (2.2) 

In [6] the authors introduced the following FALVA problem with higher-order derivatives: 
Problem 2.3. Find the stationary points of the integral functional 

nqi-)] = £ ^ ^(^), • • • , q^""^ (0)) {t - er-'d9 , (2.3) 

m > 1, under the initial conditions 

q^'Ha)=Ql, i = 0,...,m, (2.4) 

where q^{9) = q{9), q^^\9) is the derivative of q{9) of order i, T is the Euler gamma function, 
< a < 1, 9 is the intrinsic time, t is the observer time, t ^ 9, and the Lagrangian 
L : [a, b] x ]R"^(™+^) ^ M is a function of class (7^™ with respect to all its arguments. 

Remark 2.1. In the particular case where m = 1, Problem 12.31 reduces to Problem 12.11 
Theorem 12.41 generalizes Theorem 12.21 to the higher-order case. 

Theorem 2.4 ([6j). // g(-) is a stationary point of (j2.3p . then q{-) satisfy the following 
higher-order Euler-Lagrange equation.' 

d' 



^(-1)^— a,+2^(e,(Z(0),g(e),...,g(-)(e)) =F(e,g(0),(?(0),...,g(2--y , (2.5) 

i=0 

m > 1, where 



F {9, q{9),q{9), Q^'""-'^ (0)) = ^("1)^"' ^ ^^+2^ (o, q{9),q{9), • • • , 9^"^^ (^) 



k \ d^~^ 



k=2 i=2 
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We borrow from [9] the notation 



m-j 



= Y,i-^y-^d,+,+2L[9,q{e),m,---,Q^"'\0)) , j = l,...,m, (2.7) 



d9 

i=0 

which is useful for our purposes because of the following property: 



= q{9),m, • • • , 9^"^ (0)) - V^^-i , J = 1, . . . , m . (2.8) 

Remark 2.2. One can write equations (|2.5|) in the following form: 

d2L (0, qie), m, Q^'^HO)) -^i^'=F [e, q{e),q{e), q(2™-1)(0)) . (2.9) 

Theorem 2.5 ([6]). If q{-) is a solution of Prohlema \2.3\ then it satisfy the following higher- 
order DuBois- Raymond condition: 

- ^ L q{e), q{9), ■ ■ Q^""^ (0)) " E " '?^'^(^) 



= diL (0, q{e), q{9), q^""^ (0)) + F ^9, q{e), q{9), 'Z^'"'"') (0)) ■ m , (2.10) 
where F and tp^ are defined as in ()2.6p and l\2.7h respectively. 

3 Main Result 

In this work we generalize the Noether-type theorem proved in to the more general FALVA 
problem with higher-order derivatives. 

3.1 Noether's theorem for higher-order FALVA problems 

In order to generalize the Noether's theorem to Problem 12.31 (see Theorem 13.21 below) we 
use the DuBois-Reymond necessary stationary condition (I2.10|) and the following invariance 
definition. 

Definition 3.1. (Invariance of (j2.3|) ) The functional (|2.3|) is said to be invariant under the 
infinitesimal transformations 

ie = e + eT{e,q) + o{e) 
\q{9) = q{e) + eC{9,q) + o{e) 

if 

d9 



L [I m, Qio), . . . , g-'^'") (9)) it - er^'- 

= L{9,q{9)A{9),...,q^'^\9)) [t - 6) 



a-l 



+ e{t - 9r-'^ {0, q{9),q{9), • • • , 9^"""'^ (0) ) + o{e) . (3.2) 
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Remark 3.1. Expressions q'^^^ in equation p.2p . i = 1, . . . ,m, are interpreted as 

dq 



^ dd M ' ^ do^ di 

dO dO 



(i = 2, . . . ,m) . 



(3.3) 



Next theorem gives a necessary and sufficient condition for invariance of (j2.3p . Theorem 13.11 
is useful to check invariance and also to compute the infinitesimal generators r and ^. 

Theorem 3.1. (Necessary and sufficient condition for invariance of ()2.3p ) The integral func- 
tional (|2.3p zs invariant in the sense of Definition \3.1\ if and only if 

m 

diL (9, q{e),m, • • • , Q^'^HO)) T + d^+2L (o, q{e),m, • • • , Q^'^HO)) ■ 



i=0 



+ L(e,q{e),q{e),---,q'^'^Ko)]{r + ^T] = A(^9,q{e),q{e),---,q'^'"'-'Ho)) , (3.4) 



where 



dd 



{p'~^) -q^^{9)f, i = l,...,m. 



(3.5) 



Remark 3.2. If a = 1, condition I\3A\\ gives the higher-order necessary and sufficient condition 
of invariance proved in [9]: 

m 

diL {9, q{e),m, • • • , Q^'^HO)) r + ^ ^,+2^ {o, q{e),m, • • • , Q^'^HO)) ■ 

1=0 

+ L (0, q{e),m, ■ ■ 9^"^^ (0)) T = A (e, (7(0), (7(0), . . . , (0)) . 

Proof, (of Theorem 13. ip Differentiating equation ()3.2p with respect to e, then setting e = 0, 
we obtain: _ 

9 / c?*g' 



i=0 



de \de 



e=0 



+ L(r + i^r)=A. 



The intended conclusion follows from (13.3 



5e Vd0 



e=0 



d 


(d'q\ 


d 


' d 


(d'-'q\ 




de 


\de') 




de 


\d6'-^ J 


£=0- 



q^^^f, i = 2,...,m. 



□ 



Definition 3.2. (Higher-order conservation law) A quantity C {d, q{0), q{0), . . . , q^"^^ '^'> (0)) 
is said to be a conservation law if 

^c(0,q(0),(7(0),...,(?(2— i)(e)) =0 



along all the solutions q{-) of the higher-order Euler-Lagrange equation p.Sp . 



5 



Theorem 3.2. (Higher-order Noether's theorem) If the integral functional (|2.3p is invariant 
in the sense of Definition \3.f\ and T{9,q) and £,{9,q) satisfy the condition 

G {e, q{e),q{e), Q^'"^''^ (O)) ■^=-L {e, q{9),q{9), • • • , q^""^ (0)) r , (3.6) 

where 

G {e, q{0),m, • • • , Q^'""'^ (0)) = Y.{-iy-h^d,+2L [0, q{9),q{9), • • • , 9^"^ (0)) 

i=l 

+ Eg(-i)'f(S^ (, .) ^(«). m, .... (3.7) 

and n = ^ — qr, then 
C[e,q{9),q{9),---,Q^''^- 

m 

= ^ V"' + U {0, q{9) , q{9) g(-) (9)) -^^P^ ■ q^^\9) 

-K(9,q{9)A{9),...,q^^"''^\9)) (3.8) 

is a conservation law. 

Remark 3.3. Under hypothesis p.6p . the necessary and sufficient condition of invariance ([3 
takes the following form: 



diL (9, q{9),q{9), ■ ■ Q^'^HO)) r + Y, di+2L [9, q{9) 

i=0 

+ L (0, q{9),qi9), ■ ■ 9^™^ (^)) r - F (^9, q{9),q{9), ■ ■ 9^'"" 



P 



i=0 



= A(9,q{9),q{9),---,Q^'"'~'H0)) . (3.9) 
Proof, (of Theorem l3.2p We begin by writing the Noether's conservation law p.Sp in the form 

m / m \ 

C = •pO + ^V'^' L-^V'^' •g(j)(0) r-A. (3.10) 
i=2 \ i=i / 

Differentiation of equation (I3.10p with respect to 9 gives 

^ = ■ i"" + £ ('^-' i*' - i ('^-'0 

i=2 ^ ^ 

L-p^^^.qU)^e)j + lL-p^^P^.q^^\9)j . (3.11) 
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Using the Euler-Lagrange equation ()2.5p , the DuBois-Reymond condition ()2.10p , and relations 
([THI) and dTS]) in (IXTTTl . we obtain: 

m 

+ (diL + F-q)T + (^L-f^iP^ ■ q^^He)^ f 
= diLr + Lf + SaL • C + • (p^ + qr) - ■ P^ 



m 



~^^-qf + ^'^-p'^ + Y.d,+iL.pi-\ (3.12) 

Simphfication of (13.120 lead us to the necessary and sufficient condition of invariance ()3.9p . □ 

In the particular case m = 1 we obtain from our Theorem 13.21 the main result of [4j. 

Corollary 3.3. (cf. ^) If the integral functional ()2.ip is invariant under the infinitesimal 
transformations ()3.ip . and T{9,q) and^{6,q) satisfy the condition 

d3L{e,q,q)-n = -L{e,q,q)T, (3.13) 

then 

C{9, q, q) = dsL {9, q, q) ■ ^{9, q) + {L{9, q, q) - d^L {6, q, q) ■ q) r(0, q) - A {6, q, q) (3.14) 

is a conservation law (i.e., (|3.14p is constant along all the solutions q{-) of the Euler-Lagrange 
equation (12. 2|) ). 



Proof. For m = 1 we obtain from ()3.7p and (13. 8p that 

G{e,q,q)=dsL{e,q,q) (3.15) 

and 

C{e, q, q)=^^-p^ + {L{e, q, q) - ■ q{9)) r - A(0, q, q) . (3.16) 
Having in mind the equations ()2.7p and ()3.5I1 . we conclude that 



fi:'=dsL{9,q,q) , 

We obtain the intended result substituting (|3.17p into p.l6p . □ 
3.2 Example 

In order to illustrate our result, we consider an example for which the Lagrangian L do not 
depend explicitly on the intrinsic time 9. 



7 



Example 3.4. Let us consider the following second-order (m = 2) FALVA problem: to find 
a stationary function q{-) for the integral functional 

I'[q{-)] = \f^ {aq' + bq' + f) {t - Or-^dO , (3.18) 

where a and b are arbitrary constants. In this case the Euler- Lagrange equation ()2.5p reads 
5(1 -a) / (l-a)(2-a)\ / 2(1-q;)\... 

Since the Lagrangian L do not depend explicitly on the independent variable 9, the necessary 
and sufficient invariance condition (j3.9p is satisfied with 

T = l, (3.20) 
e = 0, (3.21) 



A = Fq^A = J Fqde, (3.22) 

where 

„ i — ct .... (1 — a)(2 — a) 

i. = _(4^_2,)-L_iL_Jj. (3.23) 



The conservation law ()3.8p with m = 2 takes the following form: 

C{0, q, q, q, q) = L {9, q, q, q) r + ^^sL {9, q, q, q) - -^d^L (9, q, q, ij)^ ■ n 

+ d^L {9, q, q,q)-n-A {9, q, q, q, q) . (3.24) 

Substituting the quantities L = \ {aq^ + bq^ + if), ([3:21]) and ([3:22]) into ([3H]), we 

conclude that 

i [aq^ - b(f + 3f) - q'q - J Fqd9 (3.25) 
is constant along any solution q of (13.190 . 

If a = 1, then one see from ([3.23p that -F = 0, and ([3.25p gives the classical result in [I]: 

^ (ag2 - hq^ + 3^'^) - q q 

is a conservation law. 
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